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We present a method for the introduction of small-scale structure into strings constructed from 
products of rotation matrices. We use this method to illustrate a range of possibilities for the shape 
of cusps that depends on the properties of the small-scale structure. We further argue that the 
presence of structure at cusps under most circumstances leads to the formation of loops at the size 
of the smallest scales. On the other hand we show that the gravitational waveform of a cusp remains 
generally unchanged; the primary effect of small-scale structure is to smooth out the sharp waveform 
emitted in the direction of cusp motion. 
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I. INTRODUCTION 

The study and detection of gravitational waves is a 
very active area of current research. Gravitational ra- 
diation is a direct test of dynamical gravity and eluci- 
dating its properties experimentally has the potential to 
confirm (or even change) our current gravitational the- 
ory. Furthermore, the detection of a gravitational wave 
background would open a new observational window onto 
a time in the early universe much earlier than that of 
recombination, beyond which the universe is opaque to 
electromagnetic waves. The potential rewards to cosmol- 
ogy that would result from such a detection are count- 
less. With the new interferometric gravitational wave 
telescopes (LIGO, GEO and TAMA) now operating and 
plans for a future space-based telescope (LISA) it will 
become possible to detect waveforms from various astro- 
physical sources, including topological defects. This has 
the potential to determine, or at worst constrain, the 
types of high energy particle theories that describe our 
world. 

Topological defects are a prediction of most particle 
physics models that involve a symmetry breaking phase 
transition P, as well as some brane- world scenarios Q 
and are therefore quite generic. Cosmic strings, in par- 
ticular, have drawn considerable attention because they 
do not cause cosmological disasters and are good candi- 
dates for a variety of interesting cosmological phenomena 
such as gamma ray bursts Q, gravitational wave bursts 
and ultra high energy cosmic rays P|. 

Analytic studies 0, la] and numerical simulations 9, 
ITol [Til ] show that cosmic string networks evolve toward 
an attractor solution known as the "scaling regime" . In 
this regime the energy density of the string network is a 
small constant fraction of the radiation or matter density 
and the statistical properties of the system such as the 



correlation lengths of long strings and average sizes of 
loops scale with the cosmic time t. This solution is pos- 
sible because of intercommutations that produce cosmic 
string loops which in turn decay by radiating gravitation- 
ally. 

Simulations also have found that the perturbations on 
long strings and most loops have the smallest possible 
size, the simulation resolution, which does not scale. This 
small-scale structure arises from a build-up of the sharp 
edges (kinks) produced at string intersection events 
and the initial Brownian character of the network. It is 
not clear, however, how very small loops can be formed at 
late times. The mechanism that has so far been suggested 
is that they form at regions of the string known as cusps 
|l3| , although in this case we expect the loops formed to 
have relativistic speeds. This is discussed at some length 
below. It may also be that the small loops form from the 
collapse of horizon-sized loops. The size of this small- 
scale structure is also thought to scale with the cosmic 
time t. The value typically quoted in the literature, first 
proposed in |14| . is TGfj,t, where F is a number of order 
100, G is Newton's constant and /x is the mass per unit 
length of the string, but the actual value may be many 
orders of magnitude smaller [Tsj . 

One of the key signatures of cosmic strings is the pro- 
duction of gravitational waves . The network of long 
strings and loops is expected to produce a gravitational 
wave background as well as bursts coming from cusps and 
kinks. These bursts stand out of the gravitational wave 
background and may be detectable with current sensitiv- 
ities 5] . 

In this paper we investigate the shape and gravita- 
tional waveforms of cusps with small-scale structure. In 
the next section we review the motion of strings and 
construct solutions where helical small-scale structure 
is introduced using products of rotations. In Section 
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III we investigate various small-scale structure regimes 
and show that the presence of small-scale structure at 
cusps under most circumstances leads to the production 
of loops at the size of the smallest scales. We further il- 
lustrate the various regimes, and shapes of the resulting 
cusps, using the solutions found in Section II. In Sec- 
tion IV we describe the numerical algorithm employed 
to compute gravitational waveforms and in Section V we 
show the waveforms for cusps on loops with and without 
small-scale structure. We conclude in Section VI. 



II. STRING MOTION 

When the typical length scale of a cosmic string is 
much larger than its thickness, and long-range inter- 
actions between different string segments can be ne- 
glected, the string can be accurately modeled by a one- 
dimensional object. Such an object sweeps out a two- 
dimensional surface in space-time referred to as the string 
world-sheet. Here we will consider the motion of strings 
in flat space, which is a good approximation to motion 
in a cosmological setting provided we are studying a loop 
or a segment on a long string which is much smaller than 
the horizon. 

To study the motion of strings we use the Nambu-Goto 
action |17L Il8l| which is proportional to the area of the 
string world-sheet. Minimising this action in the gauge 
where 

dux'^dux^ = dyx'^d^Xfj^ = 0. (1) 
yields the wave equation for the string 

a„a„x^(u,«) = o. (2) 

The coordinates u and v are null on the string world- 
sheet and typically given in the literature as u — t ~ a 
and V = t + a, where t is the time coordinate and a is 
a spatial parameter that measures energy density along 
the string. The constraint Eqs. leave some residual 
freedom which can be removed by imposing aP = u and 
h° = v. 

The general solution to Eq. ^ can be therefore be 
written as 

^{cT,t) = ]^[si{t-a)+h{t + a)] (3) 
and the constraint Eqs. become 

a'' = b''^l, (4) 

where a' and b' are the derivatives of a and b with re- 
spect to their arguments, t — a and t + a respectively. 
These vectors therefore live on a unit sphere 
will be referred to as the right- and left-moving halves of 
the string. 

Expressions for the functions a' and b' are most easily 
written as sums of Fourier harmonics. Unfortunately the 



unit magnitude constraint Eq. generally gives a non- 
linear set of equations involving the vector coefficients 
of the Fourier expansion and parametrising strings be- 
yond the first few harmonics proves to be analytically 
intractable. There exists a method to generate strings 
involving products of rotation matrices that greatly sim- 
plifies this task 20]. In this method the a' and b' are 
generated by starting with a unit vector and operating on 
it with a series of rotation matrices that through trigono- 
metric identities produce harmonics. The advantage of 
using this method to generate loops is that the unit mag- 
nitude constraint on the right- and left-movers is satisfied 
trivially. 

As an example of this method we can define the fol- 
lowing loop 

27r?T/ 27r7Ti 
a'{u) = Rxy{ct)Ryz{—j—u)z, h'{v) = Rxz{—j—v)z, (5) 

where Rij{^) is a 3-dimensional rotation matrix that ro- 
tates the i-j plane by an angle 7; the direction of the 
rotation is chosen such that the positive end of the ith 
axis is rotated toward the positive end of the jth axis. 
On the unit sphere, a' as given by Eq. (jsj is a unit circle 
which is tilted away from the y-z plane by an angle a. 
Similarly b' is a unit circle on the x-z plane. This loop 
is (aside from an overall spatial rotation) identical to the 
family of loops introduced by Burden (21| . 

Another advantage of the method of products of ro- 
tations is that it lends itself readily to modifications of 
the loop trajectory. It is particularly easy to introduce 
helical perturbations in one or both directions. If we ex- 
amine Eq. we observe that both a' and b' consist of 
a rotation of the unit vector z on some plane. To intro- 
duce helical perturbations in, say, b', we can replace the 
starting unit vector with 

z^z'^{l-eiy/'z + e,Rxyi^v)x, (6) 

with £{, < 1 and then subject this new unit vector to a 
rotation on the same plane as before, namely, 

h'{v) = 

„ ,27rm . 9,1/9^ „ ,27111, 

= Rxzi^vKd - e2)i/2£ + ebRxy{^v)x). 

(7) 

The unit vector given by Eq. ((HJ consists of two parts. 
The first part is a vector of magnitude (1 — e^)^^^ in the 
direction of z and it corresponds to the straight part of 
b'. The second part is a circle of radius eb on the x-y 
plane, perpendicular therefore to z, that winds around 
this plane lb times. This represents the perturbed part 
of b'. Acting on this new unit vector z' with the rotation 
Rxzi^^v) produces a curve in the unit sphere with he- 
lical perturbations of radius eb- An analogous operation 
can be carried out on a'. 

Figure n illustrates the effect of this procedure. Equa- 
tion (O gives a' and b' that trace out great circles on the 
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FIG. 1: Right and left movers on the unit sphere. The thicker 
curves that encircle the sphere once correspond to a' and b' as 
given by Eq. with a = 0, n = 1 and m = 2. The complex 
curves correspond to the perturbed versions of a' and b' with 
ea ~ ~ 0.2, la ~ 100 and h = 200. We have introduced 
twice as many oscillations in the perturbation for b' because 
the unperturbed b' winds twice around the unit sphere. 

unit sphere. Perturbations from Eq. ((2} and its analog 
for a' add small, circular oscillations to the larger loops. 



III. THE FORMS OF COSMIC STRING CUSPS 
IN THE PRESENCE OF SMALL-SCALE 
STRUCTURE 

Cusps are regions on long strings or loops that achieve 
phenomenal Lorentz boosts. These boosts, combined 
with the potentially large mass per unit length of cos- 
mic strings, may lead to the production of a detectable 
gravitational wave signal . 

Cusps arise when the right- and left-moving parts of 
the string cross on the unit sphere. At the crossing they 
point in the same direction, namely, 

a' = b'. (8) 

In this case 

x' = i[-a' + b']=0, (9) 

and 

x= i[a' + b'] =a' = b'. (10) 

Here x' = and x = dfX. The tip of the cusp therefore 
moves at the speed of light and a small region of string 



around the cusp has an enormous Lorentz boost. Cusps 
are generic in the sense that the a' and b' are confined 
to live on the unit sphere and for general string and loop 
configurations we typically expect them to cross. 

The effects on cusps of the presence of small-scale 
structure on just one of the string halves ('chiral' struc- 
ture) was investigated briefly in |22l] . There it was argued 
that the shape of the cusp depends on the contribution 
of the small perturbations to the effective mass per unit 
length of the string. Here we will review and expand that 
discussion to the more general case of small-scale struc- 
ture on both halves of the string. Although the detailed 
properties of small-scale structure in cosmic string net- 
works are still not understood it is useful to study the 
range of possibilities. Below we show that under most 
circumstances the presence of small-scale structure at a 
cusp leads to the formation of loops at the size of the 
smallest scales. 

For simplicity we consider a loop of size L with small- 
scale structure of wavelength A and amplitude to wave- 
length ratio e. For now we take the small-scale structure 
to be the same in both directions. The string has an 
effective thickness due to the small-scale structure, 

d-2eA. (11) 

When a long string or loop has a cusp the effective 
thickness of the string may result in overlap. Calcula- 
tions for the overlap of the fields making up the string 
core have been performed in 0| and 0| . In the case of 
core overlap the fields in a small section of the string near 
the tip of the cusp unwind and produce energy in the form 
of particles, the bosons that make up the string. Fortu- 
nately we can borrow these results, expecting, however, 
for overlap to result in the formation of cosmic string 
loops rather than particles. The latter work f23| took 
account of the Lorentz contraction of the core and is the 
correct calculation for generic string core overlap. The 
former work did not account for the contraction of 
the string core but, as we will see, the main result may 
be useful in the case of back-reaction-dominated cusps. 

If the small-scale structure is small enough to have 
little effect on the average motion of the string, then the 
analysis of ,'24] applies, and we expect an overlapping 
region near the cusp of length 

l^VdL (12) 

which will typically fragment into about eL/X relativis- 
tic loops of the small-scale size A. Here the effect of the 
small-scale structure has been to replace the core thick- 
ness of the string with the effective thickness d, which is 
usually much larger. 

The overlap begins when the Lorentz gamma factor of 
the cusp reaches the value 

7o ~ /LAi- (13) 

When the small-scale structure is not symmetrically dis- 
tributed between the right- and left-movers at the cusp. 
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the overlap expressions, Eqs. H12() and (|13|l . will be dom- 
inated by the larger of the two halves. 

On the other hand, if the small-scale structure is 
large enough, there will be important corrections to the 
Nambu-Goto motion. The effect of small-scale structure 
increases as the string begins to move rapidly near a cusp, 
with the effective mass per unit length becoming 

Mcff -Ai(l+7'e')- (14) 
Therefore when the gamma factor reaches the value 

7ft - 1/e (15) 

the contribution to the effective mass per unit length of 
the string becomes significant and the effective motion 
of the string deviates from regular Nambu-Goto motion. 
Averaging the motion of the string over scales of the size 
of the small-scale structure and smaller results in an in- 
crease of the effective mass density and a decrease of the 
tension of the string [25l |. 

Whether back-reaction prevents the formation of the 
relativistic cusp depends on which one of the two gamma 
factors, Eq. l(TK|l or Eq. (fT!^ . is the smallest. If e'^L/d < 
1, then 7f, > 7o, and overlap begins before back- reaction 
is important. On the other hand, if e^L/d > 1, then 
back-reaction changes the form of the cusp. 

The simplest situation of the latter sort occurs when 
both string halves contribute equally to the effective mass 
per unit length. Then the effective motion of the string 
is given by |2| 

x(ct, t) = ^ [a(fcu) + h{kv)] (16) 

with k = Vl — f^- The two halves therefore live on a 
sphere of radius fc < 1. Sharp cusps form as before, 
because at crossings on the sphere the effective x' = 0, 
but the tips of cusps move at a speed |x| = k rather 
than at the speed of light. Because the cusps are sharp 
we expect the small-scale structure to overlap. If the 
perturbations are small then the cusp is relativistic and 
we expect an overlapping region around the cusp of a size 
given by Eq. (|12|l and the production of relativistic loops. 
However, if the small-scale structure is sufficiently large 
the region of overlap will be larger due to the negligible 
Lorentz contraction of the core and we can use the result 
in p^ . In this case we expect an overlapping region of 
size 

around the cusp, which will fragment into about 
e^^^{L/ X)^/^ loops of the small-scale size A. Since the 
loops produced in this case are non-relativistic, this 
mechanism may be responsible for the large number of 
small loops seen in simulations. 

At first glance the small-scale structure in numerical 
simulations would appear too small to slow down cusps 
sufhciently to produce non-relativistic loops. However, 



the wavelength of the perturbations on long strings is 
a very small fraction of the length of the long strings 
and the fact that it can be observed at all in simula- 
tions means the amplitude to wavelength ratios cannot 
be negligibly small. 

If the small-scale structure is not symmetrically dis- 
tributed between the right- and left-movers at the cusp, 
for example because of statistical fluctuations, the radii of 
the spheres on which the effective right- and left-moving 
halves live will be different. They will therefore miss each 
other and make x' finite, with 

\^\^\k,-ka\^\el-el\ (18) 

at the cusp. Here, the sub-indices a and h distinguish 
between right- and left-movers. The limiting situation 
of 'chiral' small-scale structure (i.e., £(, = 0) was investi- 
gated in "22] and is similar to that of a superconducting 
string with a chiral neutral current ^26, , 2J\ ■ 

Since x' is finite the cusp will be rounded. If the radius 
of curvature of the cusp is comparable to the effective 
thickness of the string then we are in the same situation 
as before: a section of string will overlap near the cusp. 
The amount of overlap will be given by either one of 
Eqs. (|12|1 or (|17|) depending on the size of the small-scale 
structure. 

If the radius of curvature of the cusp is much larger 
than the thickness then there will be no overlap. This 
case is analogous to that of superconducting strings with 
a chiral neutral current and we expect self-intersections 
to occur at the cusp about half of the time . These 
self-intersections result in the production of a loop of size 

l^\el-el\'''L (19) 

which will fragment into about — e^|^/^L/A loops of 
the small-scale size A. 

To summarise, we have seen that in most cases overlap 
of the small-scale structure at a cusp leads to the forma- 
tion of loops at the smallest scales. The only exception 
occurs when the small-scale structure is a) sufficiently 
large that the averaged motion of the string deviates sig- 
nificantly from Nambu-Goto motion, and b) sufficiently 
asymmetric that the radius of curvature of the cusp is 
larger than the effective thickness of the string. Even 
in this case, however, we expect intersections to produce 
small loops about half the time. 

Figure [3 illustrates the various possibilities for the 
shape of cusps with small-scale structure described 
above. Without small-scale structure, we have a sharp 
cusp which would evaporate at the tip by core over- 
lap. With symmetrical small excitations, the cusp is still 
sharp and will chop off a loop by the overlap of the wig- 
gles. With asymmetrical small-scale structure, the cusp 
will be smoothed out. There may or may not be an in- 
tercommutation, depending on which half of the string 
has the larger perturbations, as discussed in p^ . 
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FIG. 2: On the left, a projection of an unperturbed cusp 
on a loop given by Eq. Q with q = 0, n = 1 and m = 2. 
Next to the right, a perturbed version of the same cusp with 
ta = tb = 0.3, la ~ 200 and h = 400. Next a projection of a 
cusp with Sa = 0.4, tb = 0.5, L = 200 and h = 400. In this 
case the effective x' is finite and the cusp is smoothed out. 
On the right, a cusp with n = 2, m = 1, = 0.4, et — 0.5, 
la = 200 and = 400. Switching the values of n and m has 
inverted the sign of x' and led to a self-intersection near the 
cusp. 



IV. GRAVITATIONAL WAVEFORMS: 
FORMALISM AND NUMERICAL 
IMPLEMENTATION 

Periodic sources produce a metric perturbation that in 
the wave-zone is given by a sum of plane waves p9l | 



n— — oo 

The polarisation tensor e|j"^ is given by 
4G,„ . 1 



(20) 



e(r;(x,u;„) = — (V(fc„) - ^V.uT^iK)) (21) 



and 



(22) 



is the Fourier transform of the stress-energy tensor of 
the source. In these expressions r is the distance to the 
source, T is the period of motion, 



(23) 



is the wave-vector of the gravitational wave, is a unit 
vector pointing from the source to the point of observa- 
tion and ujn = I'Kn/T is the frequency. 



In the gauge defined by Eq. |0J the stress energy tensor 
of the string can be put in the form 

T^^%x)^^l f dudv{x%x';^+x%x%)S^'Hx-xiC)). (24) 



We consider a loop of length L which oscillates peri- 
odically with period T = L/2. The Fourier transform of 
the stress-energy tensor of the loop is 

T^''(/c„) ^ j ' dt j d^xe^'^-'-^y" dudv 

X [x^„<„+<„<J<5W(x-a;(^i,^;)). (25) 

The expression for the string position in terms of right- 
and left-moving parts when substituted in Eq. H25|l yields 
the symmetrised product 

T^^'ikn) = ^{A^^{k^)B''{K) + A''(fc„)i?^(fc„)), (26) 



where 



and 



A^ikn) - Y / dua'^(M)e*'="-'^(")/2^ (27) 
^ Jo 



B^ikn) = Y f dvh'^'{v)e 
^ Jo 



ik„-b{v)/2 



(28) 



This means that we can write the polarisation tensor 
Eq. (EH) as 



(29) 



We are free to choose coordinates such that the wave- 
vector lies on the z-axis. In this case the only physical 
components of the polarisation tensor can be shown to 
be ej"' = e^' and = ^x"^ 113 1 which are referred to 
as the -I— and x -polarisations of the gravitational wave 

US-. 

It is easy to show using integration by parts and the 
periodicity of and 5^ that 



kn-^tJ.ikn) — k^^B^ikn) — 0. 



(30) 



Using this identity and Eq. (|29|) one finds that the am- 
plitudes of the -I— and x -polarised waves are given by 



(n) _ 2G'^L 



(A,(fc„)B,(fc„)- Aj,(fc„)Bj,(fc„)) (31) 



and 



^(„) ^ ^GfiL^^^^^^^^^^^^^^^ ^ Ay{k,,)B,{kn)) (32) 



respectively. 
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The metric perturbation can therefore be written as 



/i+/x(x,i)= J2 ef/^,(x,c.„)e-''="-. (33) 



The n=:0 mode in this sum corresponds to the static field 
and should not be included in the sum if we are interested 
in the radiative part of the metric perturbation. 

In order to study the gravitational waveforms produced 
by arbitrary loops we must resort to numerical methods. 
If we consider a loop composed of some continuous func- 
tions sl{u) and b(u) we can always approximate it by 
choosing some discretisation scheme that is sufficiently 
fine so that all the features of the loop are resolved. 
We can use this discretisation to compute the integrals 
Eqs. H27|l and H28I) and find the gravitational waveform 
it produces. 

Here we use the piece- wise linear discretisation |3lj . 
Allen and Ottewill 32] computed the gravitational wave- 
forms for a number of different loop trajectories and 
showed that the in the piece-wise linear approximation 
the waveform for those loops is in good agreement with 
their analytic calculations. Furthermore, as they found, 
an advantage of using this particular approximation in 
the calculation of gravitational waveforms is that the in- 
tegrals Eqs. (|27|) and H28|l can be performed exactly. This 
result is computed explicitly below. 

In the piece-wise linear discretisation both functions 
a(-u) and b(t;) are constructed out of Na and Nf, straight 
segments joined at kinks. The values of the internal pa- 
rameters u and V are continuous and are labeled by Uj 
with j — 0, 1, Na ~ 2, Na — 1, and a similar expression 
for Vj, at the kinks. Here we will take Na = Nt, = N and 
space the internal parameters evenly between the kinks 
(wj+i — Uj = Ufc+i — Vk — S for all j and k) for simplicity. 



The functions a(u) and b(u) are defined at the kinks 



aj = a{uj), bfc = b(wfe) 



(34) 



and the derivatives of these functions in the intervening 
straight segments between the kinks 



aj+i - a 



i h' - 



bfc+i - hk 




FIG. 3: Waveforms for the loop defined by Eq. (|KJ with a — 0, 
n — 1 and m = 2, whose cusp is shown in the inset and 
first in Figure |2| as a function of time. The solid line is the 
X -polarised waveform when looking straight at the cusp in 
the 2-direction. No gravity wave is emitted with + polari- 
sation in this direction. For the same loop with an overall 
spatial rotation of 1 radian about the x-axis, so that we are 
looking off-cusp, the -|— polarised waveform is dashed and the 
x-polarised waveform is dotted. All waveforms are summed 
over the first 2000 Fourier modes with 2G^L/r = 1. The loop 
is approximated with 5000 linear segments giving a resolution 
of the internal parameters 5 — 2 x 10~*I/. 



the integrand in Eq. H36|l is a simple exponential. From 
here it is straightforward to show that 



JV-l 



2a 



with a similar expression for B^{kn) 



X ^g«'^Ti(«3 + l-f^-aj + i)/2 _ giw„(uj-0-aj)/2jj-<5g-j 



(35) 



GRAVITATIONAL WAVEFORMS FROM 
CUSPS 



Here we will focus only on A^^{kn) because the situation 
for B^{kn) is completely analogous. For piece-wise linear 
loops the integral of A^(fc„) in Eq. H27|l becomes a sum 
of integrals between kinks 



L ^ J 



^^gifc,..a(«)/2^ 



(36) 



Since the string segments are are taken to be straight 
between the kinks 



a^(u) = + {u — Uj)a 



(37) 



Gravitational bursts produced by cusps have previ- 
ously been studied analytically by Damour and Vilenkin 
5]. They found the metric perturbation produced by a 
cusp to be 



h{t) (X |i-r|i/3 



(39) 



where r is the distance to the source. 

As an example we can check waveform of a cusp on a 
loop. We take the loop defined by Eq. Q which has a 
cusp in the z-direction at t = and t = L/4. The on- 
and off-cusp waveforms for this loop are shown in Figure 
|3|for a = 0, n = 1 and m = 2. 
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FIG. 4: Log-log plot of the cusp waveform for the loop de- 
fined by Eq. @ with a = 0, n = 1 and m = 2 as a function 
of time. The solid curve gives the x-polarised waveform com- 
puted numerically and the dashed curve is given by Eq. (I40II 
with a — 2.11 (fit to the data). The waveform has been 
summed over the first 2x 10^ Fourier modes with 2GfiL/r = 1. 
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(t-rJ/L 

FIG. 5: Waveforms for a perturbed loop, given by q = 0, 
n = 1, m = 2, ta = tb = 0.3, L = 200 and h = 400, 
whose cusp is shown second in Figure |5| and in the inset, 
looking straight at the cusp in the z-direction. Because of the 
perturbations there is now a small -(--polarised wave and the 
form of the x-polarised wave is not as sharp as in Figure |21 



Figure 0| shows a log-log plot of the waveform versus 
time near the cusp. The solid curve gives the wave- 
form computed numerically over the first 2 x 10^ Fourier 
modes, with 2GfiL/r = 1 and the dashed curve is given 

by 

/ix - a|t-r|i/3 (40) 

with a = 2.11 (fit to the data). The large number of 
Fourier modes necessary to reproduce the behaviour of 
Eq. (|39|l is due to the slow decrease in amplitude oi 
as a function of the Fourier mode, hx{n) oc n^^^'^. 

Due to the pervasiveness of small-scale structure in 
numerical simulations of cosmic strings, it is important to 
test the robustness of the prediction for the gravitational 
waveform of a cusp, Eq. H39|) . to the presence of small- 
scale structure. 

In order to investigate the effects of small-scale struc- 
ture on the gravitational waveforms of cusps we have per- 
formed a number of simulations of cosmic string loops 
with small perturbations. Figures 130 show the wave- 
forms of the loops whose cusps appear in Figure |21 The 
waveforms are computed by summing over the first 10000 
Fourier modes with 2G/iL/r = 1. The loops are approx- 
imated with 5000 linear segments. Intercommutations 
are not performed in the simulation; the strings just pass 
through each other without producing new loops. 

These figures show that the overall shape of the wave- 
form of the cusp is not greatly affected by the presence 
of small-scale structure even when the amplitude of the 
small-scale structure e is not particularly small. The pri- 
mary effect is to smooth out the very sharp point in the 
waveform. But that point is visible only if the observer 
is located in the direction of cusp motion. Gravity waves 
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FIG. 6: Waveforms for a perturbed loop with a = 0, n = 1, 
m — 2, ea = 0.4, e;, = 0.5, la ~ 200 and 4 — 400, whose cusp 
is shown third in Figure Inland in the inset. 



emitted from the string moving with Lorentz boost 7 are 
confined to a cone with angle 9 ~ |^. The smoothing 
of the sharp waveform is due primarily to the limitation 
of the string boost at 7{,, so the change in the waveform 
will be observed only if one is within angle 7jJ~^ '--^ e of 
the direction of cusp motion. 

In addition to the smoothing of the large cusp the 
waveforms also show significant sub-structure. It is a 
result of the small helical perturbations and is present in 
both the the -I— polarised and x-polarised waveforms. 

The sub-structure is particularly prominent at the time 
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FIG. 7: Waveforms for a perturbed loop with a = 0, n = 2, 
m = 1, £a = 0.4, €b — 0.5, Za — 200 and 4 = 400, whose cusp 
is shown last in Figure Inland in the inset. 

the large cusp forms and corresponds to the waveforms 
of many tiny cusps. These smaller cusps form as a result 
of the crossings on the unit sphere that occur around the 
time the large cusp forms (as in Figure^. In terms of 
gravitational wave detection it is possible that the small 
scale structure gives the dominant contribution to the 
waveform in some range of frequencies. This possibility 
depends on the value of G/x and the details of cosmic 
string network evolution and deserves further investiga- 
tion. 



VI. CONCLUSIONS 

The presence of small-scale structure at cusps leads 
to a multitude of interesting behaviours mostly involving 
self-intersections of the loop or long string that produces 
the cusp. The specific processes and outcomes depend 
on the detailed properties of the small-scale structure. 
Indeed, one of the scenarios presented here supports the 
ideas of Albrecht who first suggested that the small 
loops seen in simulations are produced at cusps. 



At first glance it would appear that cusps are only 
able to produce ultra-relativistic loops, which would rule 
out this mechanism since such loops are not seen in sim- 
ulations. However, in the regime where the small-scale 
structure has enough amplitude to significantly affect the 
motion of the string, the cusp and the loops it produces 
through overlap will no longer be relativistic. 

As an illustration we take the canonical value for the 
size of the small-scale structure, which for GUT strings 
IS A - IQ-H, and a significant amplitude to wave-length 
ratio, e ^ 1. Then we expect a horizon-sized cusp, of size 
L ^ i, to produce a few hundred non-relativistic loops 

ei/3(L/A)2/3 - 500. (41) 
The primary effect of small-scale structure is to smooth 
out the sharp waveform emitted by the large cusp in the 
direction of its motion. That waveform is smoothed any- 
way if one is located away from the exact direction of 
motion; observers outside a cone whose angle is propor- 
tional to the small-scale structure amplitude will see little 
change. However, since superimposed on the large cusp 
there is a multitude of smaller ones, it may be that the 
small-scale structure turns out to give the dominant con- 
tribution to the gravitational wave signal in the range of 
frequencies onto which we happen to have our observa- 
tional window. This depends on the typical size of loops 
and the spectrum of perturbations that ends up on loops. 
This possibility deserves further investigation. 

Unfortunately, the uncertainties that remain regarding 
string network evolution make it difficult to say anything 
definite in this regard. For instance, if it turns out to 
be correct that most loops are produced at the smallest 
possible scales then we do not expect them to have any 
significant structure. 
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